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Abstract. A Bose-Einstein condensate in a high-finesse ring cavity scat- 
ters the photons of a pump beam into counterpropagating cavity modes, 
populating a bi-dimensional momentum lattice. A high-finesse ring cav- 
ity with a sub-recoil linewidth allows to control the quantized atomic 
motion, selecting particular discrete momentum states and generating 
atom-photon entanglement. The semiclassical and quantum model for 
the 2D collective atomic recoil lasing (CARL) are derived and the su- 
perradiant and good-cavity regimes discussed. For pump incidence per- 
pendicular to the cavity axis, the momentum lattice is symmetrically 
populated. Conversely, for oblique pump incidence the motion along 
the two recoil directions is unbalanced and different momentum states 
can be populated on demand by tuning the pump frequency. 



1 Introduction 

Their unique coherence properties candidate Bose-Einstein condensates (BECs) as 
ideal systems to generate and probe light-atom correlations in collective light scatter- 
ing and superradiant instabilities pQ. Superradiant Rayleigh scattering experiments 
usually occur in free space |2l3j . so that scattered photons rapidly leave the inter- 
action volume, limiting the coherence time of modes propagating along the major 
condensate's axis ('end-fire modes'). On the contrary, when BECs interact with a 
high-finesse optical cavity, the correlations between scattered events can be stored in 
long-lived cavity modes [3], allowing, for instance, to study new regimes in the strong 
coupling limit [516] . Furthermore, recent experiments on collective light scattering by 
BECs in a high-finesse ring cavity have shown the possibility to employ the cavity 
sub-recoil resolution as a filter selecting particular quantized momentum states [7]. 
These experiments rely on the collective atomic recoil lasing (CARL) mechanism, en- 
visaged by Bonifacio and coworker in 1994 [8] and finally observed in Tubingen, early 
with atomic clouds as hot as several 100/iK [9] and more recently with ultra-cold 
atoms (TUJ. CARL represents the atomic analogue of the free-electron laser (FEL) 
[llj . which has been studied for a long time in Milan [12113] . In particular, Bonifacio 
and Casagrande predicted the existence of a superradiant regime in FELs |14ll5j . 
successively observed as collective light scattering in CARL [TU]. A further advance 
on CARL theory was obtained in 2001, when the semiclassical model was extended 
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to a quantum description suitable for BECs, in which the atomic motion is quan- 
tized in photon recoil momentum states |16ll7j . Next, a full quantum CARL theory 
investigated entanglement between collective momentum states and cavity modes 
[18119] . Generally, CARL is often described in a ID geometry, i.e. with pump and 
scattered modes anti-parallel and atoms recoiling after each scattering event by 2hk 
along the incident beam direction. However, in the former Superradiant Rayleigh 
scattering experiment [2] the cigar-shaped condensate major axis was set orthogonal 
to the incident laser, and two scattered beams were emitted along the condensate 
axis, with atoms recoiling at 45° with respect to the incident laser. In that case, the 
geometry was two-dimensional with two scattered end-fire modes. CARL and Super- 
radiant Rayleigh scattering in a 2D configuration have been investigated by several 
authors 20 21 22 23 24 25 26 27[2H]- More recently, a two-frequency pumping scheme 
has been implemented to enhance the resonant sequential scattering in the Superra- 
diant Rayleigh scattering [29130] and in CARL [3T]. In particular, a sub-recoil cavity 
linewidth combined with a bi-chromatic pump, with frequency separated by twice the 
recoil frequency, allows to observe subradiance in a degenerate cascade between three 
collective momentum states [32] , 

In this paper I consider a BEC in a high-finesse ring cavity scattering photons 
from a pump beam into two counter-propagating cavity modes [7] . Varying the pump 
intensity and frequency, it is possible to populate in a controlled way a 2D momentum 
lattice, where atoms belonging to different sites get entangled with the scattered 
cavity photons. The paper is organized as follow: In sec. II I derive the semiclassical 
2D CARL model for a pump beam incident at variable angle; In sec. Ill the atomic 
motion is quantized. Sec. IV discusses the different quantum regimes and present 
some numerical result in the nonlinear regime. 



2 Semiclassical model 

Let's consider N two-level atoms (with \a) and \b) upper and lower states, respec- 
tively) in a cloud with length L and diameter W <C L 7 exposed to an uniform, 
s-polarized along e y laser beam, incident in the cavity plane (x, z) and making an 
angle <j> respect to the normal of the cavity's optical axis z (see figHJ, with electric 
field 

JfJ. _ 2l. ^J^^ e i(k x cos <p-k z sin <p-ui t) _|_ c c | (J) 

with ujq = cko- The pump photons are scattered in two counter-propagating cavity 
modes with frequency coinciding with a cavity eigenfrequency, uj c — ck c , and electric 
field 

E s = y jj^eKfeeS-We*) + E2(} -i(k e Z+ Uc t) + c c | 

The pump and the two cavity mode fields induce the following coherence between 
the states \a) and \b), 

Pab = \ {s< oe l Mz cos 0- Z sin 0-ct) + g^ik^z-ct) + g^k^z+ct) + Q Q j (3) 

and a force F = d y V(Ei + E s ) y in the cavity plane (x, z), where d y = d(p a b + c.c.) 
is the y-component of the electric dipole moment and d is the dipole matrix element. 
Assuming the pump-atom detuning A a = luq — cu a much larger than the spontaneous 
decay rate r, it is possible to see that Si — S7i/A a (for i = 0, 1,2), where ]?,; = 
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Fig. 1. 2D CARL configuration: the pump field Eo is incident with an angle (ft with respect 
to the normal of the cavity axis z, and the cavity modes Ei and E2 are counter-propagating 
along z. All the fields are linearly polarized perpendicularly to the cavity plane. The atoms 
recoil along the ei and directions, when they scatter Eo pump photons into the modes 
Ei or E2, respectively. 



dEi/h 8 . A straightforward calculation shows that the equations for the momentum 
components p x z = mv x z are [23] 



dp x 



dt 4A a 
dp z . hko fio 
~~dt ~~ ~ l AA a 



cos0|^ie" iqix + r? 2 e" 



c.c 



} 



sin<£ |i7ie" jqi x - J? 2 e~ iq2X - c.c.} 



2A„ 



} 



where we assumed k c ~ fco, we introduced fi\^ 2 = ^1,2 exp(iA c t), where A c = ujq — uj c 
is the pump-cavity detuning, and 

qi,2 = fc [cos^e x - (sin^± l)e z ], (6) 

where qi • q 2 = 0. The equations for the cavity mode amplitudes are 

dfl\ ckod 2 n a 



dt 2ie{)TiA a 
dQ 2 ckod 2 n a 
dt 2it^hA a 



{f2 Q (e lcllX ) + + A>(e^ 2lfe ° 2 )} - k«A +*AA (7) 
{/2 (e tq2 ' x > + ^2 + /)i(e 2lfe ° 2 )} - k c Q 2 + iA c f2 2 (8) 



where n a is the atomic density and k c = cT/L c is the linewidth of the ring cavity 
with length L c and transmission T . It is more convenient to describe the atomic 
motion along the directions of qi^ = 91,261,2 with unitary vectors ei 2 and gi,2 = 
k§\j2{\ ± sin (f>). Then, the momentum components along these directions are, in 
units of the photon recoil momentum hqi^ 2 , 



Pi, 2 = — = ^-3— [Pxcos0 -p z [sva(p± l)\ 



Hqf 



Defining the phases 



9i,2 = qi,2 • x = kox cos (ft — &oz(sin 4> ± 1) 



(9) 
(10) 
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and the dimensionless field amplitudes 



where V is the interaction volume, the complete equations for N atoms and the two 
cavity mode amplitudes are 

— = 2u rlPlj 

de 2j 



dt 



2ul r2 p 2j 

aia 2 e 

^.g {02e -.«.. + , c . } _,^ {ai(l|e .<.„-» ul _. 





(12) 




(13) 




(14) 




(15) 




(16) 


" aJ" 2 


(17) 



where j = 1, . . . ,N, g = \J d 2 loq / {2e KV) is the single-photon Rabi frequency and 

Wri,2 = i (1 ± sin0)cj r , (18) 

where w r = 2?ifc 2 /m is the maximum photon recoil frequency. Notice that the atoms 
move along the recoiling directions e\ or e 2 (see figfT]) when they scatter the pump 
photons into the cavity modes a\ and a 2 , respectively (see the first terms on the right 
hand sides of Eas. (fl"4|) and (fl"5|) . representing the dipole forces depending on 9\j and 
9 2 j, respectively). Furthermore, the atoms recoil further along the cavity axis z when 
they exchange photons between the two cavity modes themselves (see the second 
terms on the right hand sides of Eas. (fl4| and (fl"5j). representing the dipole force due 
to the two-cavity mode interference and depending on 8 2 j — 9±j — 2koZj). Notice 
that the longitudinal dipole force breaks the pump-atom detuning A a symmetry: In 
fact, if A a — > — A a and — > — dip, these terms change sign too (together with the 
collective single-photon light shift Ng 2 /A a ). 



3 Quantum model 

In a quantum theory, the classical variables 9±j , 9 2 j , pij , p 2 j , «i and a 2 are promote 
to operators, with commutation rules [9 a j,P0m] = iS a p5j m and [o a , ot] = 8 a p where 
a,/3 = 1,2. Without cavity losses (i.e. n c = 0), Eas. l|T2"j) - l[T7|l derive by the following 
Hamiltonian: 



A' 



H = 1Z 



3 = 1 



2 i 2 , -9^0 

u} rl p l3 + UJ r2 p 2j + i^- 



g 

A a 



( A , t t X 

I A c - I (a[ai + a 2 a 2 ). 



(19) 
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The single-particle Hamiltonian Hi (Oi , 62 , p\ , P2 , Q>i , cl\ , 0.2 , a\ ) can be second-quantized 



MS 



p2rr p2-rr 

H= Mi de 2 ^\e 1 ,92)H 1 {e 1 ,92,-idg 1 ,-id 62 ,a 1 ,a\,a 2 ,al)W{e 1 ,9 2 ) (20) 
Jo Jo 

where the quantum field operator <?(#i,02) obeys bosonic equal-time commutation 
rules [&(9 1 ,9 2 ),&(6[,6 2 )} = <5(0i - 9[)8{9 2 - 9' 2 ) and f{9 1 ,0 2 ),4f(9 l l ,9 l 2 )] = 0, with 
normalization condition J dB\ J d9 2 ^'(9\, 02)^(^1, 62) — N. Introducing the anni- 
hilation operators for the two momentum components pi and P2, i.e. ^(^1,^2) = 
Y^ m ,n^rn, n u m (0!)u n (9 2 ), where u m (0i, 2 ) = (1/V^tt) cxp[im6»i, 2 ] and [c m , n , & m , „,] = 
8 m ,m,'8 n ,n' , the Heisenberg equations for c m! „, d\ and 0,2 read: 



dc r , 



dt 



dt 



-i [m 2 u} r i + n 2 Lo r2 ] c m . n 



2A a 
.9 2 



d&i _ gftp \ - 



Cm— 1 ,n+l 



dd 2 _ gf2p \ - „ t . _ -t 

_ 2Z\ y C m,n c ™,"-1 * a a l / j c m 

a m,n a m,n 



(21) 
(22) 

(23) 



where A — A c ~ Ng 2 /A a . Notice that in Eq. (j!?Tj) we have neglected the global phase 
factor proportional to A(a\a\ + a\a,2)- 

In the following we will neglect the quantum nature of the operators c mj „ and 
ai,2 and we treat them as complex dynamical variables. Furthermore, we introduce 
dimensionless time, r = 2uj r pt, field amplitudes, A\^ — (ai^/VpN) exp(-iAt), and 
pump parameter, Aq = a^/^/pN, where eto = S7o/(2g) is such that Oq is the pump 
photon number. The collective CARL parameter p is defined as [5] 



apgWN 
2A a w r 



2/3 



(24) 



Then, defining C m ^ n — (l/viV)c m , n exp[i(m + n)At], k = n c /(2u r p) and S = A/u) r , 
Eqs.{|ll])-(|23j) yield' 



dC-rn^n 

dr 



dAi 
dr 



I 

2p~ 



1 + sin < 



1 — sin ( 



— (m + n)5 



+ p [A\C m -i. n + A*2C m ,n-\ — AiCm+i.n — ^CWn+l] 

— i~7- [AiAjCm+i^^i + A* l A 2 C m -l,n+l\ 

Ao 

— J]] C* n n C m -i^ n — 1— C*j . n Cm-\,n+\ — K,A\ 
m,n ^ m,n 

Ax 



An 



(25) 
(26) 

(27) 
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Notice that the total probability is conserved, i.e. 'Yl, mn \G v 



1. The growth 



rates for the two modes are Gi,2 = — 2Im(A)(2o; r p), where A is the solution of the 
cubic dispersion relation: 



-5 



Sl,2 

2p 



+ 51,2 = 0, 



(28) 



where si 2 = (l±sin0)/2. Notice that the longitudinal lattice term is nonlinear (being 
proportional to Ai^) and it does not contribute to the dispersion relation ([28]) . In 
the following we will indicate the 2D momentum lattice states as (to, 71), associated 
with momentum components pi = m(hqi)ei and P2 = 771(^,(72)62, respectively. In 
the linear regime, the two cavity modes grow independently and atoms (initially in 
(0,0)) populate the states (±1,0) and (0, ±1), respectively. In particular, when the 
laser beam is parallel to the cavity axis z (i.e. tj> = 90°), s\ = 1 and s 2 = 0: Only the 
mode a\ grows, so we can set 0,2 = and the model reduces to the usual ID CARL, 
with c m , n = S nfi c m [IS]. 



4 Discussion 



After scattered a photon with momentum Kk and energy fuj, the atom recoils with 
momentum p determined by energy and momentum conservation laws, i.e. ftko = 
Kk + p and Hujq = hui + p 2 /2m, so that p x = hkg cos <fi, p z = hko^simp =p 1) and 
lu = ojq — (j r (l ± sin0)/2, where the upper and lower sign is for a photon emitted 
in the cavity mode ai or 02, respectively. Since u> is near the cavity mode frequency 
u c , then by tuning the pump frequency u>o near the cavity mode frequency ui c it is 
possible to enhance one mode with respect to the other, depending on incidence angle 
<fi, gain bandwidth ui r p and cavity linewidth k c values [7]. 



4.1 'Good-Cavity' and 'Superradiant' regimes 

The cubic dispersion relation (|28|) provides the expression for the gain rates G1.2 of 
the two cavity modes in the good-cavity (GC) regime (i.e. when 2n c -C G\^) and 
in the super-radiant (SR) regime (i.e. when 2k c > Gi^), either in the semiclassical 
regime (i.e. when Au> 3> u> r ) or in the quantum regime (i.e. when Auj < oj r ), where 
Auj is the resonant gain bandwidth [16117118] . In particular, in the quantum regime 
the atoms scatter the pump scattered photons only forward, since the atomic recoil 
red-shifts the scattered photon frequency (i.e. at w - ^,1,2) such that it is set out 
of the resonant gain bandwidth 16J. So, in the quantum regime the atoms populate 
initially only the positive-momentum states (1, 0) and (0, 1). In the GC limit we can 

set k c — 0, obtaining Gi,2 = 2^J (2p) 3 uj 2 — A\ 2 for A\£ < (2p) 3 / 2 w r , where = 
A — w r i,2 is the detuning taking into account the recoil shift. Hence, in the quantum 
GC regime, the maximum gain and the gain bandwidth are G max = 2(2p) 3 / 2 uj r and 
Auj — 2{2pf'l 2 uj r — G max , respectively, and the conditions necessary to observe it 
are p < 1 and k c <C uj t . In the quantum SR regime, Gi^ = K c (2u) r p) 2 / (Af 2 + k 2 ), 
so that the maximum gain is G max = (2u> r p) 2 / k c and the gain bandwidth is n c [17] . 
The condition necessary to observe the quantum SR regime is ui r p < k c < u> r . Notice 
that in the quantum regime gain and bandwidth are the same for the two modes. 
However, increasing p for a given cavity linewidth n c , the system moves toward the 
classical GC limit, Gi,2 ui r) where the recoil shift can be neglected and the gain is 
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Fig. 2. Gain coefficients g± t 2 = —2Im(X) (continuous blue line for mode 1 and dashed 
red line for mode 2) vs. the pump-cavity detuning 5 — A/u) r , for <j> = 45°, p = 0.4 and 
k c = 0.8oi r . 



centered around A = 0, with G max — 2\/3 w r p(si ) 2) 1 ' 3 - Hence, in the classical regime 
the two cavity modes have different gain rates. As an example of an intermediate 
case (with parameters close to those of ref.[7]), figH] shows the gain G\ t 2 (in CARL 
bandwidth 2pui r units) vs. the pump cavity detuning S — A/uj r , for <f> — 45°, p = 0.4 
and k c = 0.8tJ r . 



4.2 Symmetric nonlinear regime 

The case where the laser beam is perpendicular to the cavity axis z (i.e. <f> = 
and 51,2 = 1/2) has been discussed in details in ref.[23]. Here, the two modes are 
symmetric and atoms move at 45° forward and backward with respect to the cavity 
axis. This is also the original configuration of the Superradiant Rayleigh scattering 
experiment of ref.[2]. Following ref . |24j . it results that in the quantum regime atoms 
populate sequentially the momentum states (n,n) (with n = 1,2,...) by a four- 
level 'diamond' transition, passing through the intermediate states [n — l,n) and 
(n, n — 1). Since c„ jm = c TOi „ and A\ — A2, each transition (n — 1, n — 1) — » (n, n) 
can be described by optical Bloch equations for two-level systems once a popula- 
tion difference W — |c n _i.„_i| 2 — |c„ jn | 2 and a polarization S = c*_ x n _ 1 c n _i ) „ + 
c n-i n c n.n are introduced. In the quantum SR regime and at resonance (i.e. for 

A = cj r /2), the populations evolve as |c„.„(t)| 2 = (1/4) |l — tanh[-^/ p/k(t — ro)]| 

and |c„_i !rl (r)| 2 = (l/4)sech 2 [-\/ p/k(t — td)], whereas the cavity photon number is 
(n p h) = iVsech 2 p/ k(t — td)], where td ~ yfp\n(p). Notice that the maximum 
occupation probability of the intermediate states (n— 1, n) and (n, n— 1) is 1 /4. This 
configuration is particular attractive, since either entanglement |33j and subradiance 
could be there easily addressed, as discussed by Crubellier et al. [34135] . 



4.3 Asymmetric nonlinear regime 

The symmetry of the case with perpendicular incidence is broken when the laser 
beam shines the atoms with an oblique incidence angle, as can be seen for instance 
in figlU In this case, changing the laser frequency ujo with respect to the cavity mode 
frequency u> c allows to unbalance the two counter-propagating cavity modes, as well as 
the two momentum components pi and p2- Furthermore, nonlinearity induces more 
complicated dynamical structures resulting from the interplay of cooperative gain 
and cavity losses, when more than one photon is scattered by the condensate. As an 
example, we consider the case with <f) — 45°, p — 0.4, k — 0.8w r and different pump- 
cavity detuning A. In order to ge the analysis simpler, we neglect the longitudinal 
lattice (i.e. the last term in the right-hand side of Eas. (l25]) and the second terms in 
the right-hand side of Eqs. (|2"o]) and (|57|) assuming Ai, 2 <C A Q . Fig[3]shows the result 
of numerical integration of Eas. (l25[) - (|27p for different pump-cavity detuning values: 
A = —2.5oj r (left column), A = — 0.025a; r (central column) and A = 2u) r (right 
column); mode intensities |^4i,2| 2 and the average momentum components (pi^) are 
shown vs. r in the upper and lower lines, respectively; blue thick lines refer to mode 
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XXX 

Fig. 3. |^4_i,2 1 2 (upper line) and (pi,2) (lower line) vs. dimensionless time r for A — — 2.5u r 
(left column), A = — 0.025u r (central column) and A = 2uj t (right column), for (j> — 45°, 
p = 0.4 and k c — 0.8lj t - Blue thick lines are for the mode 1, red thin lines for the mode 2. 



1 and red thin lines refer to mode 2. For A = — 2.5w r (left column) G\ <C G%, so 
that the mode 1 does not grow appreciably. The atoms move along the 62 direction, 
up to the momentum state (0, 2) after a time r = 180. For A = — 0.025o; r (central 
column) the gain rates are equal (G\ = G2 = 0.41w r ) and the atoms initially equally 
populate the momentum states (1,0) and (0,1). However, later on the atoms turn 
toward the 62 direction, populating the states (1, m) with m — 1, ... ,6 after a time 
r = 180. This rather peculiar behavior has been observed also in the experiment 
of ref.[7]. This behavior can be easily understood observing that the recoil shift for 
the mode 1 is s\/s2 = 5.83 times larger than for the mode 2, so that the incident 
photons are set out of resonance after the atoms have scattered the first laser photon 
into the mode 1, whereas the incident photons remain well inside the resonant gain 
bandwidth when scattered into the mode 2 (see figEJ). As a consequence, the atoms 
scatter a single photon into the mode 1, stopping at the momentum states (m, n) with 
m = 1, whereas they are allowed to scatter photons into the mode 2 (up to n = 6, 
as shown in central column of fig'E]), populating the momentum states (l,n). Finally, 
for A — 2oj r (right column of figEJ) G2 <C G\, so initially only the mode 1 grows 
and atoms populate sequentially the states (1,0) and (2, 0); however, at a longer time 
the mode 2 grows and reaches saturation, so that the atoms populate also the state 
(2,1). 



5 Conclusions 

I have derived the semiclassical and quantum model of the collective atomic recoil laser 
(CARL) for a Bose-Einstein condensate set in an arm of an high-finesse ring cavity, 
with a laser beam incident at an oblique angle with respect to the cavity axis. The 
atoms scatter photons into two counterpropagating cavity modes, recoiling along two 
different directions determined by the incidence angle. For perpendicular incidence, 
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atoms scatter symmetrically the pump photons into the two cavity modes, populating 
sequentially symmetric momentum states (n,n) with n = 1,2,.... Conversely, for 
oblique incidence it is possible to populate in a controlled way different momentum 
states by tuning the pump frequency, as experimentally done in ref.[7]. Similarly to 
the ID geometry |18ll9j . it is expected that atoms belonging to different momentum 
states may be entangled between themselves and/or with the photons scattered in the 
cavity modes. Furthermore, an even reacher scenario can be realized by a bichromatic 
pump with frequency spacing tuned around the recoil frequency, which is expected to 
enhance or inhibit the transfer of atoms between different momentum states [31132) . 
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